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Abstract 

We study cosmological braneworld models with a single timelike extra di- 
mension. Such models admit the intriguing possibility that a contracting 
braneworld experiences a natural bounce without ever reaching a singular 
state. This feature persists in the case of anisotropic braneworlds under some 
additional and not very restrictive assumptions. Generalizing our study to 
braneworld models containing an induced brane curvature term, we find that 
a FRW-type singularity is once again absent if the bulk extra dimension is 
timelike. In this case, the universe either has a non-singular origin or com- 
mences its expansion from a quasi-singular state during which both the Hubble 
parameter and the energy density and pressure remain finite while the cur- 
vature tensor diverges. The non-singular and quasi-singular behaviour which 
we have discovered differs both qualitatively and quantitatively from what is 
usually observed in braneworld models with spacelike extra dimensions and 
could have interesting cosmological implications. 
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Following the seminal papers [1,2], most braneworld models with extra dimensions as- 
sume that extra dimensions are spacelike, so that the brane is embedded in a Lorentzian 
multidimensional manifold. However, there is no a priori reason why extra dimensions 
cannot be timelike, and the observational constraints on such braneworld models were dis- 
cussed in [3]. More recently, such models were under consideration in [4-8]. In this letter, we 
demonstrate that a timelike extra dimension could, in the case of the simplest braneworld 
model, lead to interesting new features. For instance, a contracting braneworld generically 
bounces as it reaches a high density thereby leading to the absence of the 'big crunch' and 
'big bang' singularities of general relativity. We consider both Friedmann-Robertson-Walker 
(FRW) and Bianchi I scenarios and demonstrate that, under some additional and not very 
restrictive assumptions, the presence of anisotropy does not modify this conclusion. 

In this letter, we consider the case where a four-dimensional hypersurface (brane) is the 
boundary of a five-dimensional Riemannian manifold (bulk) with nondegenerate Lorentzian 
induced metric. The action of the theory has the natural general form 

S = M 3 \ [ (K - 2A) - 2e / k] + [ (m 2 R - 2a) + / L(h ab , 0) . (1) 

. J bulk J brane J J brane ^ ' J brane 

Here, 1Z is the scalar curvature of the five-dimensional metric g a b in the bulk, and R is 
the scalar curvature of the induced metric h a b = g a b — en a nb on the brane, where n a is 
the vector field of the inner unit normal to the brane. The quantity K = K ab h ab is the 
trace of the symmetric tensor of extrinsic curvature K ab = h c a V ' c n b of the brane. The 
parameter e = 1 if the signature of the bulk space is Lorentzian, so that the extra dimension 
is spacelike, and e = — 1 if the signature is (— , — , +, +, +), so that the extra dimension is 
timelike. Since the induced metric on the brane is assumed to be Lorentzian, the signature 
of the extra dimension coincides with the type of the unit normal n a to the brane. The 
symbol L(h a b, 0) denotes the Lagrangian density of the four- dimensional matter fields <fi the 
dynamics of which is restricted to the brane so that they interact only with the induced 
metric h^. All integrations over the bulk and brane are taken with the natural volume 
elements \J—eg d 5 x and \J — h d 4 x, respectively, where g and h are the determinants of the 
matrices of components of the metric in the bulk and of the induced metric on the brane, 
respectively, in the coordinate basis. The symbols M and m denote, respectively, the five- 
dimensional and four-dimensional Planck masses, A is the five- dimensional cosmological 
constant, and a is the brane tension. In this paper, we use the notation and conventions of 
[9]. 

Variation of action (1) gives the equation of motion in the five-dimensional bulk: 

Gab + Agab = , (2) 

and on the brane: 

m 2 G a b + crh ah = eM 3 S ab + r ah , (3) 

where Gab an d G a b are the Einstein's tensors of the corresponding spaces, S a t, = K a ^ — Kh a b, 
and Tab denotes the four-dimensional stress-energy tensor of matter on the brane. 
Using the Gauss relation on the brane, one obtains the constraint equation 

e (R - 2A) + K ab K ab - K 2 = . (4) 
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Then, expressing the extrinsic curvature K ab from (3), one obtains the following scalar 
equation on the brane [7]: 

eM 6 (R - 2A) + (m 2 G ab + ah ab - r ab ) (m 2 G ah + ah ab - r ab ) - ± (m 2 R - 4a + rf = , 

(5) 

where r = h ab r ab . 

Before proceeding further, we pose to discuss the significance of the timelike character 
of the extra dimension. We emphasize that, in the theory under consideration, the extra 
dimension is accessible only for the gravitational field, while all matter fields are constrained 
to the brane and are propagating on the background of the induced Lorentzian metric. Thus, 
the dynamics of the matter fields on the brane has the standard properties of quantum 
field theory in curved Lorentzian spacetime, and the effect of the extra timelike dimension 
tells only in the gravitational sector. An important issue demanding further investigation 
is related to the tachyonic nature of the Kaluza-Klein gravitational modes that could in 
principle lead to violation of causality and unitarity on the brane. Some discussion of this 
issue within the context of braneworld models with more than one time like dimension can 
be found in [3,4]. 

In passing, we mention that a generalization of the Randall-Sundrum solution [2] to the 
case of arbitrary signature of the extra dimension and arbitrary Ricci-flat vacuum brane can 
be written in the form [4,7] 

ds 2 = edy 2 + exp [~^^v) h a p(x)dx a dx p , (6) 

where h a p(x) are the components of the Ricci-flat metric on the brane, which is situated at 
y = 0, and the bulk coordinate is in the range y > 0. Equation (6) is supplemented by the 
constraint A e ff = 0, where the expression for A e fj is given by Eq. (10) below. One can see 
that, for e = — 1 and negative brane tension a, the gravity appears to be 'localized' at the 
brane (it is, perhaps, more appropriate to speak of 'deflation' in the case of timelike extra 
dimension) . 

In this letter, we study cosmological implications of this braneworld theory. We first 
consider a homogeneous and isotropic cosmological model with the cosmological time t, 
scale factor a(t), energy density p{t), and pressure pit). In this case, by integrating Eq. (5), 
one obtains the following equation [7]: 



0, 



where C is the integration constant, H = a/a, and k = 0, ±1 corresponds to the spatial 
curvature of the universe. For m — 0, which corresponds to the Randall-Sundrum limit [2], 
this equation reduces to 



H 2 , « = A , ^ 2 , 2arp ep 2 C 

a 2 6 9M 6 M 6 M 6 a 4 ' { ' 



Introducing the notation 
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A ea 2 . . 

A =»=2 + 3M^ (10) 




one can write Eq. (8) as follows: 

rr2 , K A e ff , 8irG N p ep 2 C 

Thus, G N is the effective gravitational constant, and A eff is the effective cosmological constant 
on the brane. (For e = 1, Eq. (11) reduces to the well-known results [10,11].) 
A FRW universe described by (11) is embedded in the bulk with the metric 

fdv 2 

ds 2 = - f(r)du 2 + ^- + r 2 dn K , (12) 
f( r ) 

where 

/ A„2 /nr\ 

(13) 

and dQ K is the standard Euclidean three-dimensional metric corresponding to k, which 
satisfies the vacuum bulk equation (2). The position of the brane is described by the equation 
r = a{u). 

The function /(r) is assumed to be positive in the domain of action of the coordinate r, 
which, in particular, implies that at least one of the constants k, A or C must be nonzero. 
Metric (12) has singularities and horizons in the bulk for certain values of the constants. We 
will not concern ourselves with this issue, assuming that such singularities can be avoided, for 
instance, by the introduction of another brane. Although the presence of another brane can 
affect the evolution and spectrum of perturbations (as in the case in the Randall-Sundrum 
model), it does not modify the general equations of the cosmological evolution of the brane. 

We first discuss some interesting consequences of Eq. (8) or (11) for e = —1. Specifically, 
as we briefly noted in [8], the fact that the p 2 term on the right-hand side of (11) has a 
negative sign, leads a contracting universe to bounce instead of reaching the singular state 
p — > oo, R a bcdR abcd — > oo — typical of general relativistic 'big crunch' singularities. In order 
that the subsequent evolution of the universe be compatible with observations, one requires 
the brane tension a to be negative so that the effective gravitational constant, given by (9), 
is positive. The additional requirement that the bounce take place at densities greater than 
that during cosmological nucleosynthesis leads to the constraint [11] |er| > (lMeV) 4 . Since 
e = —1, a braneworld with A e ff > in (10), must have a positive bulk cosmological constant 
A > 0. 

As the universe collapses to high densities, the negative p 2 term grows much faster than 
the remaining terms on the right-hand side of (11) leading to H — > and to the inevitable 
bounce of the braneworld. This feature appears to be quite generic, since it depends neither 
upon the equation of state of matter nor upon the spatial curvature of the universe. The 
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simplest singularity-free bouncing braneworld model (with C, k, and A e g equal to zero) has 
the form 

3 M 6 y ! 

An example of a bouncing universe containing radiation is shown in Fig. 1. It is easy to 
show that a spatially closed universe (k — 1) with matter satisfying p+3p > will be 'cyclic' 
in the sense that it will pass through an infinite number of nonsingular expanding-collapsing 
epochs. As demonstrated in [12], a massive scalar field in such a universe usually leads to 
an increase in the amplitude of consecutive expansion cycles and to a gradual amelioration 
of the flatness problem. 




FIG. 1. A bouncing radiation-dominated braneworld (rj = J dt/a is the conformal time). 

In passing, we note that cosmological braneworld equations with p 2 correction terms of 
negative sign on the right-hand side were also recently discussed in [13] in the context of 
the Randall-Sundrum two-brane model with a bulk scalar field stabilizing the radion. We 
also note that bouncing and cyclic braneworlds were discussed in [14] in the context of the 
usual theory with spacelike extra dimension where the bulk metric is that of a charged anti 
de Sitter black hole. In this case, the bounce occurs if the black-hole charge is sufficiently 
high. In our model, the bounce is quite a generic feature. 

Next, let us consider the case of a homogeneous but anisotropic braneworld described by 
the Bianchi I metric 

3 „ 

ds 2 = -dt 2 + J2a 2 i (t)(dx i ) . (15) 
i=i 

In this case, Eq. (5) with m = gives the following closed equation on the brane: 

6H + 12H 2 + £ (Hi - Hf = 4A cfr + [a(p - 3p) - p(p + 3p)] , (16) 
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where we used the notation of (10) and 



H = ^, H=Wh i = ~, a=( ai a 2 a 3 ) 1/3 . (17) 
The last term on the left-hand side of Eq. (16) is the shear scalar of the spatial section 

3 2 

of the universe: {Hi — H) = cr a p(J a 1 . Because of the presence of this term, and because 

the evolution of the shear tensor cr a p is not specified on the brane (see [15]), Eq. (16) cannot 
be integrated completely, but the result of its integration can be written in the form 

H-^jtw-HUM^ + ^ + g + Z, (18) 

1=1 

where the notation is the same as in (9) and (10). We emphasize that Eq. (18) describing 
the evolution of the Bianchi I brane is absolutely general. 

As we noted, the behaviour of the second term on the left-hand side is not specified 
and, in principle, can be arbitrary. Nevertheless, our conclusion about the bounce of the 
contracting universe in the case of a timelike extra dimension (e = — 1) remains valid as long 

3 

as the shear scalar cr a/3 a al3 = J2 (Hi — H) does not grow faster than a~ 8 as a — > during 

i 1 

the contraction of a radiation dominated universe. As before, the bounce is caused by the 
negative p 2 term on the right-hand side. However, in the case of the anisotropic model, the 
bounce is taking place only as regards the overall expansion a/ a, while the behaviour of 
the shear tensor remains unspecified in general; this leaves open the possibility that such a 
universe bounces only in a given spatial patch. 

As an example of bouncing behaviour, the shear can be specified [15,16] by the additional 
assumption that the so-called nonlocal energy density on the brane U, which is determined 
by the projection of the Weyl tensor to the brane, either vanishes or is negligible. This leads 
to [15] 

3 «V2 

a a/3 a^ = ]T (Hi - Hf = — , £ = const , (19) 
i=\ a 

which clearly reinforces our earlier conclusion about the bounce. Indeed, the integral in 
Eq. (18) can now be evaluated to give 

2 A cfr SvrGV ep 2 C E 2 

H = — + ^r + w« + ^ + ^> (20) 

and it is immediately clear that the p 2 term on the right-hand side, which grows as a" 8 during 
the radiation-dominated contraction stage, dominates over the last (shear) term. The same 
result is obtained in a slightly more general setup [17], if the nonlocal energy density U is 
assumed to behave like radiation, U oc a -4 . Note that, although the universe bounces, the 
behaviour of the individual components of the shear tensor a af3 (hence also of the individual 
components Hi, i — 1, 2, 3) is not fully specified in the above examples. 

The collapse of a scalar-field dominated universe is even more likely to lead to a bounce 
because the energy density p^ of the field <f> during contraction becomes kinetic-dominated, 
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2 / 2 » V(<j>), so that p<f, oc a 6 , which makes the p 2 term in (11) grow much faster than in 
the case of radiation [12]. 

So far, we have been assuming m = in the action (1). Let us return to a homogeneous 
and isotropic universe but drop this constraint, i.e., extend our study to braneworld models 
whose action contains the induced curvature term on the brane. The cosmological equation 
(7) corresponding to this model can be solved with respect to the Hubble parameter: 



H 2 + 



K 
^2 



P+ (J 

3m 2 



2e 

72 



1 ± Wl + e£ 2 



p + a 
3m 2 



A 
6" 
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A C e 

- + -T + 



I 2 



1± \ l + e£ 2 



p + a 
3m 2 



A 
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C_ 
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where we introduce the length parameter 



(21) 



I 



2m 2 
IP 



(22) 



Similar equations were obtained in [18] for the case of e = 1. The 'i'signs in (21) correspond 
to two different ways of bounding the bulk space by the brane, depending on whether the 
inner normal to the brane points in the direction of increasing or decreasing bulk coordinate 
r in (12). Alternatively, the two different signs in (21) could correspond to the two possible 
signs of the five- dimensional Planck mass M. The model with the '— ' sign is the one that 
passes smoothly to the previously considered case of m = in the sense that the right-hand 
side of (21) with '— ' sign can formally be expanded in powers of m. The right-hand side of 
(21) with '+' sign formally diverges as m — > 0. 

Since model (21) with '— ' sign passes smoothly to the model described by Eq. (11) as 
m — > 0, it is clear that, for sufficiently small values of m, the evolution described by (21) 
with '— ' sign will not be very different from that described by the limiting Eq. (11) with the 
same initial values of the scale factor and energy density, and the bounce will take place in 
the case of e = —1. This is explained by the fact that the last term in the second line of (21) 
has a negative sign and grows by absolute value during the contracting phase. In particular, 
the bounce will definitely take place in the case of spatially flat or closed universe with zero 
dark-energy term (C = 0) if M 2 /6 < 1 =>• Am 4 /M 6 < 3/2. If the value of m is not so 
small as to lead to the bounce of the contracting universe, then such a universe will reach a 
singularity similar to that of the case e = 1 described in our paper [19]. Specifically, as the 
universe collapses, its energy density p grows and can, under certain assumptions, exceed 
the dark radiation term under the square root in (21). Further increase in p will cause the 
expression under the square root to reach zero, heralding the formation of a cosmological 
singularity beyond which the solution cannot be continued. This singularity is unusual, since 
the energy density and pressure remain finite, while the space-time curvature of the brane 
and its extrinsic curvature tend to infinity as the singularity is approached (see [19] for more 
details). 

Our general conclusions are therefore the following: 

1. The four-dimensional cosmological evolution of a Randall-Sundrum-type model with 
a timelike extra dimension is non-singular both in the past and in the future provided matter 
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satisfies the weak energy condition p + p > 0. The 'big crunch/big bang' singularities are 
also absent in a broad class of anisotropic Bianchi I braneworld models with a timelike extra 
dimension. 

2. The presence of an 'induced' brane curvature term in the higher-dimensional action (1) 
can trigger the formation of an initial or final singularity even when the extra dimension is 
timelike. However, the properties of this singularity are unusual since the Hubble parameter 
remains finite while the curvature tensor diverges as the singularity is approached. 

Thus, braneworlds with timelike extra dimensions have properties which are fundamen- 
tally different both from standard general-relativistic behaviour and from the properties 
of braneworld models with spacelike extra dimensions, and which could have interesting 
cosmological implications. The important issue of the tachyonic nature of Kaluza-Klein 
gravitational modes in such theories and the related problem of unitarity has been discussed 
in [3,4] but demands more extensive examination. 
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